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Abstract 

We discuss several issues related to a recent proposal for denning classical spatial av- 
erages to be used in the so-called cosmological backreaction problem. In the large 
averaging-volume limit all gauge dependence disappears and different averages can be 
univocally characterized by the observers associated with different scalar fields. The 
relation between such averaging procedure and the standard one is emphasized and 
a gauge invariant way to select different observers is presented. For finite averaging 
volumes we show that, within our proposal, a residual gauge dependence is left, but is 
suppressed by several effects. 



1 Introduction 



The geometric properties of the Universe, on sufficiently large scale of distances, may be 
described by a homogeneous and isotropic Friedman-Lemaitre-Robertson- Walker (FLRW) 
metric. On the other hand on small scales the Universe exhibits a very inhomogeneous 
structure. In this way the FLRW metric should be interpreted as an "averaged" cosmo- 
logical metric. A problem (see, for example, [1]) appears due to the fact that the Einstein 
equations for the averaged geometry are different, in general, from the averaged Einstein 
equations. In fact the averaging procedure does not commute, in general, with the non- 
linear differential operators appearing in the Einstein equations. As a consequence, the 
dynamics of the averaged geometry is affected by so-called "backreaction" terms, originat- 
ing from the dynamic contribution of the inhomogeneities present in the metric and matter 
sectors. 

Following the discovery of cosmic acceleration on large scales, interest in these issues has 
considerably risen (see [2] for a review). In fact, it has been suggested (see e.g. [3]) that the 
dynamical effects of the backreaction could replace dark-energy sources in the explanation 
of such cosmic acceleration, solving, in this way, also the well-known "coincidence problem" . 

In order to have a physical interpretation the results that one obtains in the investi- 
gation of the backreaction problem should be independent from the gauge chosen. More 
generally, the gauge issue consist in extracting physically meaningful results from cosmo- 
logical perturbation theory eliminating all possible gauge artifacts. At linear order in the 
perturbations, the way to do this is well known: following Bardeen's pioneering work [4], 
it consists of defining gauge invariant combinations of the perturbations themselves. The 
problem is more involved when one is interested in evaluating backreaction effects either 
at the classical or at the quantum level. Such a problem has been faced in [5], where a 
gauge invariant prescription to average in cosmological setup was defined, and in [6] where, 
using the previous prescription, a general-covariant and gauge invariant generalization of 
the effective equations (introduced in [7]) which describe the cosmological backreaction was 
given (see [8] for a specific application of the prescription introduced in [5]). The key point 
of our approach is to define a general formula for the classical or quantum average of any 
scalar quantity, on hypersurfaces on which another given scalar quantity is constant. The 
physical object is the hypersurface defined by the above scalar, and not the gauge chosen 
to calculate such an average. 

The gauge issue, in connection with the backreaction, has been already partially dis- 
cussed in the literature (see e.g. [9]-[ll]), in particular in [11] the authors, after introducing 
a generalization of the effective equations presented in [7] , applied the effective equations 
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to two different gauges. Following [5, 6] the results of [11] should be seen in the following 
way. As already noted in [5] the standard average prescription present in the literature [7] , 
applied in a particular gauge, can be seen as the gauge invariant prescription introduced 
in [5, 6] with the average performed on the hypersurface defined by a scalar field which 
is homogeneous in that particular gauge. So, from this point of view, it is clear that dif- 
ferent gauges should lead to different results for the standard averaging procedure, since 
this correspond to averages performed on different hyper surf aces. The different results ob- 
tained in [11] actually correspond to the cosmological backreaction with respect to different 
observers/hypersurfaces. A key question comes out: which hypersurface should be chosen 
to perform the average? The point is to understand what are the physical observers, with 
respect to whom the backreaction should be calculated for the problem under consideration. 

In this paper, after establishing a correspondence between the standard average pre- 
scription and the one introduced in [5, 6], we will see how to define, for any given gauge, 
an associated homogeneous scalar. We will then show how the physical properties of the 
observers associated with such scalars are actually independent from the gauge chosen to 
perform the calculation. In particular we will present a gauge invariant description of ob- 
servers in geodesic motion. Such observers have been often suggested as the physically 
relevant ones (see, for example, [9]) for the calculation of the cosmological backreaction of 
the present space-time inhomogeneities. Furthermore we will discuss in detail the resid- 
ual gauge dependence of the average prescription introduced in [5, 6], when the average is 
performed on a finite volume (see Sect. 2 of [5]). 

The paper is organized as follows. In Sect. 2 we recall the averaging prescription intro- 
duced in [5, 6]. In Sect. 3 we give some useful relations among the scalars used to define the 
hypersurfaces where the average is performed. In particular the expression which defines 
the scalar that is homogeneous in a given gauge, and the relation which connects different 
homogeneous scalars in different gauges. Furthermore, we define the geodesic motion of an 
observer, up to second order in perturbation theory, in a perturbed FLRW space-time, and 
we introduce the hypersurface which defines such observers in geodesic motion indepen- 
dently from the gauge considered. In Sect. 4 we study the residual gauge dependence of the 
averaging prescription as a function of the size of the volume of integration. Our conclusive 
remarks are presented in Sect. 5. 

2 Averaging prescription: a recent proposal 

To evaluate the dynamical contribution of the inhomogeneities, present in the metric and 
matter sectors, one has first to introduce a well defined averaging procedure with respect to 
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the observer involved in the measurements. For this purpose let us begin from the following 
four dimensional integral of a scalar S(x) defined in [5] 

F(S,n)= [ d 4 xJ-g(x) S(x) = f d 4 xJ-g(x)S(x)W Q (x). (2.1) 
Jn(x) v Jma 

The integration region Q C .M4 is defined in terms of a suitable scalar window function 
Wn, selecting a region with temporal boundaries determined by the space-like hypersur- 
faces T,(A) over which a scalar field A(x) takes constant values, and with spatial boundary 
determined by the coordinate condition B < tq, where B is a (positive) function of the co- 
ordinates with space-like gradient d^B, and tq is a positive constant. As we are interested 
in the variation of the volume averages along the flow lines normal to £(A), we choose the 
following window function [6]: 

W n (x) = u"V^(A(x) - A )e(r - B{x)) (2.2) 

where 9 is the Heaviside step function, and = —jzqFJ^t^ — yj2- 

Using the window function (2.2) we can define the averaging of a scalar S(x) as 

F(S,n) = Jd 4 xV^Su»V^(A(x) - A )6(ro - B(x)) 
\ Mo.ro p^fy Jd 4 x^=g~u»V tl 6(A(x)-A )e(r -B(x))' l ' j 

In order to explicitly rewrite the above definition in terms of spatial (three-dimensional) 
integrals we can perform the covariant derivative and exploit the properties of the conse- 
quent delta function. On considering the coordinate transformation to the bar coordinates, 
in which the proper time t goes to i defined by t = h(i, x) and where the function h is 
chosen to make the scalar field A homogeneous, i.e. 

A(h(t,x),x) = A(i,x) = A(°\t) , (2.4) 

one can perform the time integration and obtain 

js Ao <PxyJ\*j{t ,x)\ S(t ,x)8(r - B(h(t ,x),x) 
{S)Ao ' ro = j^T d^xVm^Wl 0(r - B(h(t , x), x)) 

where we have called to the time i when A^°\t) takes the constant values Aq and we 
are averaging on a section of the three-dimensional hypersurface Sa„, hypersurface where 
A{x) = Aq. In this notation, for example, S is the variable S transformed to the coordinate 
frame in which A(x) is homogeneous. In this framework the hypersurface A{x) = Aq is the 
one which identifies our observers (see Sect. 3). 

As discussed in [5], if B is not a scalar, the spatial boundary can be a source of breaking 
of covariance and gauge invariance. In the limit of large averaging volume with respect 
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to the typical scale of the inhomogeneities, all the gauge dependence disappears ([5]) and 
the classical average, or the correspondent vacuum expectation value (see [5]), is uniquely 
characterized by the hypersurface of integration A[x) = Aq regardless from the gauge chosen 
to perform the calculation. In Sect. 3 we will give some useful tool to define the scalar 
A(x) and, using a particular example, we will see as the physical properties of the observers 
associated to such scalar are independent from the gauge chosen, while in Sect. 4 we will 
discuss in more quantitative terms such a breaking. 



3 Gauges vs observers in a gauge invariant averaging pre- 
scription 

As stated in [5, 6] and in the introduction, backreaction effects are usually calculated, in 
the literature, in different gauges using the standard averaging prescription [7]. In this way 
the results can be seen as particular cases of the average prescription recalled in Sect. 2 only 
if the gauge considered is the one where the scalar A(x), which defines the hypersurface on 
which we perform the average, is homogeneous. In such a gauge, using the ADM formalism, 
the standard averaging prescription [7] is also covariantly defined. In [7] this corresponds to 
identify A(x) with the scalar potential of the fluid and choose the so-called covariant fluid 
gauge [12]. 

To better understand the connection above let us see in detail how it is always possible 
to define a particular scalar which is homogeneous in a particular given gauge. 



3.1 Scalars/Gauges correspondence 



We shall work in the context of a spatially flat, FLRW background geometry, and we expand 
our background fields {g^u} up to second order in the non-homogeneous perturbations, 
without fixing any gauge, as follows: 



5oo = - 
9ij = a' 
+ 



1 _ 2 a « - 2«( 2 ), 



S i3 (l - 2^) - 2^)) + D^EU + E^) + \ ( X g + xff + *g>) 

1 ( Y (2) + y (2) , h m 



(3.1) 



where Djj = didj — 5ij(V 2 /3), and a = a(t) is the scale factor of the homogeneous FLRW 



metric. Here V> (1) , are pure scalar first-order perturbations, and 



X 



(i) 



are transverse vectors (d t B^ = and d l x\ 1> = 0), h\j' is a traceless and transverse ten- 



,(!) 



(1) 
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sor (flifcg) = = h\ 1]i ), and the same notation applies to the case of the second-order 
perturbations. 

From Eq.(3.1) one obtains 10 degrees of freedom (which must be added to those coming 
from the matter sector) which are in part redundant. To obtain a set of well defined 
equations (Einstein equations + equations of motion of the matter sector), order by order, 
we have, for example, to set to zero two scalar perturbations and one vector perturbation. 

The choice of such variables is called a "gauge" choice and can be performed using 
the gauge transformation, associated to a "infinitesimal" coordinate transformation, from 
the general gauge (3.1) to the particular gauge chosen. Such "infinitesimal" coordinate 
transformation is parametrized by the first-order, e^, and second-order, e^, generators as 
[13] 

x^x» = x» + e f 1} + \ (efoft,^) + ef 2) ) + . . . (3.2) 

where 

£ (1) = ( £ (l)' e 7l)) ' e (2) = ( e (2)' g (2)) ( 3 - 3 ) 

and we can explicitly separate the scalar and the pure transverse vector part of , e^ 2 ) as 

g (i) = ^ e (i) + e (i) ' ^(2) = ^ e (2) + e (2) • ( 3 - 4 ) 

Under the associated gauge transformation (or local field reparametrization) - where, by 
definition, old and new fields are evaluated at the same space-time point x - a general tensor 
changes, to first and second order, as 

T (i) f{i) = T« - L e(1) T(°) (3.5) 

T (2) ^ f(2) = T (2) _ Le(i)T W + I (l 2 (i) T(°) - L e(2) T(°)) (3.6) 

where L €(1 2) is the Lie derivative respect the vector 2 ^ . In particular the scalar metric 
perturbations change, to first order, as 

a (1) =a (1) -^), (3.7) 

^ {1) =/3 (1) -^? 1) + 2ae (1) , (3.8) 

^ (1) = ^ 1) +ffej> 1) + ±V 2 e (1) , (3.9) 
E<V = E w - 2e (1) . (3.10) 
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and to second order, as 



a 



( 2 ) = a^- 



where 



2 



+2 { e m e (i),i 



a 
~2 



^> = ^) + £ e 2) + V £(2) - e » 1) (2iJ«« 1 > + «« 1 >)-f^ 1) 



if 

1 1„ 3 d l &> 



(i) 



In* 
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^ = fi(2) -^-2^iw ni3 ' 



(3.11) 
(3.12) 

(3.13) 
(3.14) 



(i) 



?(!) 



1 



-«( 



(i) ^ 2 « 



+ 2^ (1) 



ae (i)i e (i) 



-2ae(i) 



D 



(i) 



(3.15) 



n l7 = -2e?^Ai 



(i)^i 
o 



+ 



V 2 

+2He° (l) (e {1)iJ + e (1)j)i ) + 



^ e (l),i e (l),i 



2^ (1) (ecijij + e^ij 
+e(i)fc,iefi)j + 2 e ?i) ( e "(i)iJfc + " 2 g ( fe i) (xSjk + *S* + ^Si) 

1 „fc - -fe \ 1 ( (1) (1) - k 

2 (, e (i)i,fce(i),i + e (i)i,fc e (i)jJ - g (Xi,k + Xfc,i + ^ J e (i)j 



+ 



+ 



For the vector metric perturbations one obtains, to first order 



M 1} 



+ 2ae (1)l 



_ (i) _ 2 

Aj — Aj Zt (l)« 



(3.17) 
(3.18) 
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and to second order [14] 



o(2) = o(2) + 



-(2) (2) 

Xi = Xi - e (2)i 



Of 



(2)i 



(3.19) 



(3.20) 



The metric tensor perturbation is gauge invariant to first order (h- = h\- ) while to second 
order one obtains [14] 



VI 



VI 



did j , 



VI 



V 2 +0; 



VI 



QkQl 



2 

V 2 " 



ii,, ■ ^ (did k u jk + d 3 d k u ik 



(3.21) 

Let us note that the differences between Eqs. (3.9,3. 10)-(3. 13,3. 14) and corresponding ones 
of [14] are connected with a slightly different definition of the scalar perturbations associated 
with gij. 

To conclude the associated gauge transformation of a scalar field A(x) = A^\t) + 
A^\x) + A^ 2 \x) is, to first order, 



A W ^ A^=A^-e° {1) A^, 

and, to second order, 

A m ^ i(2) = A (2) _ e o i)i( D _ ^ + d i e ^ d . A (i) 



(3.22) 



1 

+ 2 



(3.23) 



In this way the general scalar v4(x) homogeneus in the gauge chosen (identified by the 
generators ej^ and ej^) will be given by 

A<®(t) = A(°)(t) 

we have a class of scalars defined by their homogeneous value A^°\t) which can take any 
value. From another point of view, the coordinate transformation, that goes from a general 
gauge to the bar gauge where A(x) is homogeneous, fixes the corresponding transformation 
parameters as follows: 



= (i) i(o) ' 



; (2) 



A& am AW 



i ^(^) +e «)^ (1) ' (3 - 25) 



7 



as it is easy to control the two set of equations (3.24) and (3.25) are one the inverse of the 
other. 

To check the results in Eqs. (3.24) one can consider the simple case of a space-time 
sourced by a single scalar field 4>(x) = (f>(°\t) + (p^\x) + <p( 2 \x), and choose the uniform 
field gauge (UFG), namely the gauge for which such scalar <j>[x) is homogeneous, as the 
gauge where our scalar A(x) should be homogeneous. In such a case is pretty obvious that 
A(x) is given by a general function of 4>{x), namely A = A(4>(x)). If we now expand up to 
second order in perturbation theory we obtain 

v; m " <j><PV 0(0)^ 2^°) 2 \A(°) <f>(P) V ' 

which is, in fact, the expected solution of Eqs. (3.24). 

As pointed out the scalar A{x) is used to identify our observers. They are the ones 
sitting on the spacelike hypersurfaces over which the scalar field A{x) takes constant values, 
and their motion is given by the comoving 4-velocity 



d^A 



{-d a Ad°Afl 2 

which takes the general form 



(3.27) 



4°) = (-1,0 



a 



2 e (i),i e (i) - 2 e (i) e (i).* " a e (i)^ ~ ~r~ " 2 ra/wj + J (3 - 28) 

This normal vector to the hypersurface is the one which describes the physical proper- 
ties of our observers. As we will see in a particular example in the next subsection, such 
properties are the same in all possible gauges, the vector describes a gauge invariant set of 
properties. Let us also put in evidence that such a vector is independent from the homoge- 
neous value ^4(°) (t) of the corresponding scalar. In general, this value can be chosen at will 
and does not influence the physical properties of the observers and the physical meaning of 
our averages. For example, we can choose A^°\t) = t to have standard derivatives in the 
gauge invariant formulation of the cosmological backreaction [6] . 

An interesting result is that all these particular scalars, homogeneous in different gauges, 
can be connected to each other by a simple formula. Let us call with Ai and Aj the general 
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expressions, without fixing any gauge, of the scalars homogeneous in the gauge i and j, one 
obtains the following result 

± = A, - (Af) t - (Af) i - (e^ id M?)i (3-29) 

where (A^)i is the value that A^ takes after a gauge transformation from the general 
gauge to the gauge i (so it is gauge invariant for construction), and ( e (i)) is the value that 
we have to give to the parameter to go from a general gauge to the particular gauge i. 
It is easy to understand that this last quantity changes in the following way under a gauge 
transformation identified by a generator 

( e (i)) l -( e (i)) l - e (i) < 3 - 3 °) 
so the quantity obtained by using Eq.(3.29) is homogeneous in the gauge i and changes 
as a scalar up to second order. Starting with a scalar homogeneous in a given gauge we 
can easily obtain a scalar homogeneous in any other given gauge using Eq.(3.29). This is a 
general formula that can be generalized to any order. 

To conclude we have that any given gauge has a class of scalar fields, homogeneous in 
that gauge, while for any given scalar there is a class of bar gauges (the parameters of the 
coordinate transformation which gives the bar gauge are only partially determined by A(x), 
see Eq.(3.25)). 



3.2 A privileged observer? 

One of the central problems in the calculation of backreaction effects from inhomogeneities is 
to choose the hypersurface with respect to which our observers are defined. Such a choice will 
obviously depend on the particular backreaction problem considered. In the particular case 
of the backreaction of large scale structure on the present cosmological evolution different 
authors (see, for example, [9]) have claimed that such physical observers are given by the 
observers in geodesic motion. This property should be a gauge independent property of 
our observers, i.e. the results of the backreaction should depend only on the hypersurface 
which defines the observers and not on the gauge chosen to perform the calculation. 

An observer in geodesic motion is described by a velocity vector (normalized as 
v^v^ = — 1) which satisfies the following equation 

g aT v a V TVfl = 0, (3.31) 

using the perturbed FLRW metric given in Eq.(3.1) we obtain the following result up to 
second order in perturbation theory 

„w> = (-1,0) 



9 



= U 2 \ (3.32) 



where 



-i + Bf ) (fl (1W + BW') - a< 2 > (3.33) 

Let us note that this vector will take the standard zero order form (—1,0) in the so- 
called synchronous gauge (where = = 0, (3^ = and = 0). This suggests that 
the scalar which defines the "right" hypersurface should be the one which is homogeneous 
in such a synchronous gauge. Such scalar field A(x) can be constructed starting from 
Eqs.(3.24) or using the relation in (3.29) and the scalar homogeneus in the UFG given in 
(3.26). After some calculation we obtain that this is given by 

A (o) {t) = A (o)(t) 

A^(x) = l -A^^Jdta^ 2 + A^ j 'dt[-±aW 2 + ±(J ^a?) (/ ' dta^ 

~h (/ dta ^) + B(1)i ) + 1 + B H + B(1)i ) + a(2) ] (3 - 34) 

and we can easily check that 

u^ = v^ , (3.35) 

namely the observers defined by this scalar are always in geodesic motion independently 
from the gauge chosen. This is a gauge independent property namely a physical property 
of the observers. 



4 Gauge dependence for finite averaging volume 

In the calculation of the backreaction effects the region of integration has a typical size that 
goes from tens of Mpc to the Hubble radius. Consequently, it is important to investigate 
how the average introduced in [5, 6] depends on the gauge chosen as a consequence of the 
finiteness of the region of integration. The gauge dependence, up to second order in the 
cosmological perturbations, of a 4-dimensional integral, using the window function (2.2) 
and the coordinate trasformation (3.2), is given by (see, also, [5]) 

F(S, n) - F(S, il) = - [ d 4 x^ S(x)u^V ^(Aix) - A ) 5 (r - B{x)) 
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dB(x) ^ l dB(x) 
dx» e (!) 2 dx» 



e (l)^ e (l) ' e (2) 



(4.1) 



where a ~ indicates gauge trasformed quantities. Using this result in the averaging prescrip- 
tion defined in (2.3) one obtains: 



( s )a , 



ro 



dB(x) „ ldB(x 
dxf* e W + 2~dx^ 



F(S, n) - J d A x^) Six^V ^6(A{x) - A )5 (r - B(x)) 



1 d 2 B(x) „ „ 
2dx»dx» e ^ e ^ 



(i) Uut (i) "•" fc (2) 
F(l, 0) - | d 4 x^ n^V M ^(^(x) - ^ )<5 (r - B(x)) 



gg(g) „ 1 gg(a) 
<9a^ + 2~dxy 



V r\ (X I M 



1^B(^ 

2 5^WW 



-l 



(4.2) 



Keeping only terms up to second order in the cosmological perturbations, one obtains, after 
some algebra, the following result for the gauge dependence: 



J Saq d 3 x S {1) e(r - B(t , x)) d*x Jgg^ 5(r - B(t , x)) 
J Sao (Px 0(r o - B(t , x)) J Sao cPx 0(r o - fl(t , x)) 

k Ao d^ S^^^jro-B^x)) 



(S) 



Ao,ro 



fx Ao (Px6(r -B(to,x*)) 



(4.3) 



where is the first order vector generator of the gauge transformation of a bar quantity. 
Namely, if we consider a scalar and vector quantity in the bar coordinate and we apply 
a gauge transformation we obtain new quantities that can be obtained by the general 
definition of a gauge transformation with ej^ substituted by ej^. In practice we obtain 
that 

efi) = (0,%)). (4.4) 

Let us remark some points about this result: first, the result is independent from the window 
function used to carry out the classical average. If we repeat the procedure starting from 
Wq(x) = S(A(x) - A )6(r - B(x)) instead of u^V fl O(A(x) - A )6(r - B(x)) we obtain the 
same result. In general, we obtain the result in Eq.(4.3) starting from a window function 
5(A(x) — Ao)6(ro — B(x))C{x) with C{x) any scalar field with a non-vanishing homogeneous 
value. Second, the gauge dependence goes to zero for ro — > +oo. Third, regardless of the 
scalar S the gauge dependence is always at least of second order, so this will be always 
subleading if in the average {S)A ,r there is a non zero first order contribution. And, in 
the end, if S« = the average is gauge invariant up to second order. 
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4.1 An example 



Let us now consider a simple example where we take a sphere as spatial boundary. Namely 
we take B(to,x) = (x 2 + y 2 + z 2 ) 1 ^ 2 = \x\. Having set the spatial boundary our result will 
depend only on the structure of and e^. In general, if we consider a cosmological 
background sourced by an energy-momentum tensor with no transverse part for if , namely 
with no term Tf such that d % Tf = 0, we have that the first order vector perturbations can 
be set to zero [15] regardless of the gauge and we can set, from Eq.(4.4), = (0,d t e^). 

Using these assumptions Eq.(4.3) can be simplified to 



{S)A ,r - (S)a 



f d 3 x S W 6(r - \x\) J d 3 x V\x\ • Ve (1) 5{r - \x\) 
fd*x S W V\x\-Ve (l) 5(r -\x\) 

4^ ( 4 - 5 ) 

where we have neglected the suffix Y±a ■ If we now change coordinate system and we move 
to the spherical coordinates (r, 6, <p) where 

V/(r, 0, 0) = f-e r + -%e e + (4.6) 
or r od r sin ocp 

we obtain the following result 

~ / # d6 dr r 2 sin S {1) 6{r - r) J dxj> d6 dr r 2 sin 9 ^e {l) 5(r Q - r) 

W)A ,ro ~ W)A ,ro ~ 4^3 4^3 

/ d<j> d9 dr r 2 sin 6 S {1) |-e (1) S(r - r) 



(4.7) 



If ^ and we take a radius of integration ro which is small compared with the typical 
size of our cosmological perturbations r* ~ 100 Mpc, the integral of such first order per- 
turbation is non zero and, as said, the gauge dependence will be subleading with respect to 
this first order contribution. 

Let us now consider a region of integration with ro > r* and make the hypothesis that 
we can set to zero the integral of first order quantities. Then, following [16], we connect our 
fluctuations with Gaussian primordial ones with zero mean and we use ensemble average 
([17, 18]) to calculate the expectation value of our averages. Following the standard notation 
we denote this additional average by an over-bar. 

Let us consider the momentum expansion of a general first order perturbation <tW in 
the form 

= (^372 / d 3 ke^ s a^E(k) . (4.8) 
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Using such expansion and the properties mentioned above in the Eq.(4.7) the first order 
integrals disappear and one obtains 

-i 



(S)A ,r — {S)A ,r — 



1*4 



1 



(27T) 



d(j) d9 dr r sin 6 



I 



d 6 k 1 d 6 k 2 e 



iki-r 



o(l) 



S%E(h) ik 2 cos n e ik2 -% )k2 E(k 2 ) 5(r - r) 



(4.9) 



where Q is the angle between k 2 and r. If our fluctuations are statistically homogeneous, 
then E is a unit random variable satisfying E*(k) = E(—k) and 



E{k 1 )E{k 2 ) = 5{k l + k 2 
It follows that the ensemble average of Eq.(4.9) gives 

'4 



(4.10) 



((S) Ao - (S)a ) 



-^4 



COS £1 € 



(l)k 



(27T) 

S(r 



dO dr r 2 sin 8 



s (1) 

d k -g(&) lk 



(4.11) 



where O is the angular separation between the vector k and f. By considering the spatial 
integral inside the momentum integral, it is easy to check, by symmetry arguments, that 
the spatial integral vanishes and the gauge dependence goes to zero as a consequence of 
the assumptions made. This is not surprising, the assumption that our average over a 
limited region may be replaced with an ensemble average is reasonable only if we perform 
the integration on a region which is much bigger than the scale of our inhomogeneities. In 
this case the gauge dependence should, indeed, vanish 1 . 



5 Conclusion 

This paper is focused on different issues associated with the evaluation of the backreaction 
effects of the classical (quantum) inhomogeneities, as the residual gauge dependence of the 
averaging prescription introduced in [5, 6] and the choice of the observers with respect to 
whom we evaluate the effects. 

In Sect. 3 we have described how to construct a scalar field A(x), homogeneous in a 
given gauge, which defines the hypersurface on which we want to perform the average 
thus introducing the tools needed to have a correspondence between the standard average 
prescription [7] and the one introduced in [5, 6]. Such a scalar will identify the observers 

1 Actually, if one applies an ensemble average directly to Eq.(4.7), using Eqs.(4.8) and (4.10) and without 
setting to zero the integral of first order quantities, the first term on the r.h.s. can give a non-vanishing 
contribution which is however negligible in the limit r S> r» . 
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with respect to whom the average is performed. The physical properties of such observers 
are independent from the zero mode of the scalar field and from the possible choice of gauge. 
We have a gauge invariant way to select the observers. As a key example we gave the scalar 
and the velocity vector which identify the observers in geodesic motion with respect to the 
perturbed FLRW space-time. 

The results obtained in Sect. 3 can also be used for the evaluation of the quantum 
cosmological backreaction. In fact, the classical averages performed over all 3-dimensional 
space can be replaced by quantum expectation values using the correspondence illustrated 
in [5]. 

In Sect. 4 we have considered and formalized the residual gauge dependence associated 
with the averaging prescription (introduced in [5, 6]) when the average is performed on a 
finite volume. Let us recall again the main results: first, the result (4.3) is independent 
from the particular window function used. Second, we confirm (as already stated in [5]) 
that the gauge dependence goes to zero for averages over all the 3-dimensional space. Third, 
regardless of the averaged scalar S, the gauge dependence is always at least of the second 
order. And, in the end, if = the average is gauge invariant up to second order. 
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